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We have studied the dependence of the electronic energy spin-splitting of 
Ga^Ini-zAs-Alylni^As based double quantum wells (narrow gap structures) un- 
1 der in-plane magnetic fields. To do this, we have developed an improved version 

of the Transfer Matrix Approach that consider contributions from the abrupt inter- 
faces and external electric fields when tunneling through central barrier exists. We 
have included the dependence of the Lande (7-factor on the external applied field. 
Variations of the electron spin-splitting energy lead to some peculiarities of the den- 
sity of states. Because the density of states is directly related to photoluminescence 
excitation, these peculiarities can be observed by this technique. 

PACS numbers: 72.25.-b, 73.21.-b 
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I. INTRODUCTION 

In the last decade a great interest has arisen for the so-called spintronic, or spin-based 
electronics. The reason is that the spintronic, due to its low power consumption, promises 
to be a good alternative to the traditional electronics, based on the charge transport 1 . On 
the other hand, the spin transport is not very dissipative, with low energy losses over long 
distances, although the relaxation of the spin polarization may exist. Spintronics is a catchall 
term that refers to the potential use of the spin rather than the charge in electronic devices. 
For it to be useful in practice, spin up and spin down electronic states (or hole states) of any 
material must be separated in energy. Also, the material should be electrically polarized as 
in conventional electronics, which means that carriers, both negative (electrons) and positive 
(hole), must be able to conduct. However, except giant magnetoresistance (GMR) and its 
use in spin valves, spintronics is still just a promising possibility. 

A key point is the choice of the material. For example, narrow gap semiconductors 
with strong spin-orbit coupling. GaAs-based materials seem to be the suitable candidates 2 
due to the long life of the magnetic spin state of photoexcited electrons, which behaves 
coherently. It is also important that the materials have a good lattice matching to avoid 
defects accumulation at the interface and internal strains, which would worsen the transport 
of the polarized spin. There exist techniques to relax this strain, as to place step graded 
buffers with a progressive variation of the In concentration between the substrate and the 
active device. These layers mainly absorb strain and defects. 

Another important point is a large Lande factor for having a significant splitting of spin 
states by applying small magnetic fields. In general, In x Gai- x As — In y Al\_ y As structure 
seems to be one of the most appropriate for spintronic purposes. This heterostructure offers 
the possibility of manipulate the gap width, the Lande factor and the interface contributions 
by varying x and y concentrations. Besides, it presents a remarkable spin-splitting energy 
when a weak magnetic fields is applied. This peculiarity makes this material suitable for 
high temperature spin-valves devices. 

In this work we will focus in the spin-splitting changes in asymmetric double coupled 
quantum wells ACQW as function of an in-plane magnetic field. For this purpose we will 
base on an extended version of the 6x6 Kane formalism to calculate the band structure 
and dispersion laws^. Such electronic dispersion laws (and the corresponding spin-orbit 
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splitting) in quantum heterostructures have been widely studied in the last decades™ In 
theory, effects from compositional parameters and external perturbations can be analyzed 
through the density of states. And the modifications of the density of states can be directly 
observed using photoluminescence excitation technique (PLE) 5 . 

In bulk materials spin-orbit interaction is caused either by a soft potential 6 and by cubic 7 
and linear- spin-dependent contributions to the effective Hamiltonian. However, in the 
two-dimensional (2D) case, we can reduce the cubic contribution to a linear one after the 
squared momentum substitution by its quantized value due to confinement 9 . Moreover, it 
is necessary to consider the additional spin-orbit splitting caused by the interaction with 
abrupt heterojunction potentials (see Ref .iSLlQili) . This contribution is absolutely different 
to the contributions mentioned above. Splitting of 2D spin states is not feasible without 
the aforementioned short-range potential contribution. Unfortunately, experimental studies 
have not thrown enough light on induced mechanisms affecting the relative contributions of 
bulk and heterojunctions 6 - - -. 

The effect of an in-plane magnetic field on the energy spectrum in nonsymmetric het- 
erostructures is caused by the Pauli contribution to the electron Hamiltonian. Several pecu- 
liarities for transport phenomena in heterostructures have been also discussed^— . However, 
these discussions only consider the mixing between Pauli contribution and effective 2D spin- 
orbit interaction. As mentioned above, we cannot forget the contribution of the abrupt 
barriers in asymmetric structures, as ACQW can be, to the spin polarization. The present 
calculations, based on the three-band Kane model with nonsymmetric boundary conditions, 
try to contribute with some light to these interactions. Besides, we have included possible 
changes of the Lande factor. Although the Lande g-factor depends on the applied fields*^ 7 -, 
this dependence is negligible for in-plane magnetic fields and low electron density. The 
model also applies to narrow-gap heterostructures whenever the slow potential generated 
by doping or external transverse electric fields can be described self-consistently. 

II. EIGENSTATE PROBLEM 

If we keep in mind the strain effect (in the case that there exists) is included in the 
structure through the gap, the conduction and the valence well potentials 18 , we will center 
our attention in the Hamiltonian describing the electronic behavior in the conduction band. 
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We will include effects of electric and magnetic fields as well as the interfaces contribution. 

Based on the assumption that the gap energy e g is smaller than the energy distance be- 
tween the valence band (f-band) and the spin-split band extrema, the electronic states in 
these narrow-gap heterostructures can be described by the three-band Kane matrix Hamil- 
tonian 

e(z) + (v-7r), tt = p-^A, (1) 

where the kinetic momentum, tt, contains the vector potential A = (Hz, 0,0), H||Oy is an 
in-plane magnetic field and p = (p,p z ) is written in the p, ^-representation through the 2D 
momentum p. We have also introduced the diagonal energy matrix e(z) whose elements fix 
the positions of the band extrema and the interband velocity matrix v. 

From now on it is necessary to introduce a new index // = w, b to denote narrow-gap 
regions (wells) and wide-gap regions (central barrier and lateral sides), respectively. In the 
parabolic approximation we can write the Schrodinger equation for ACQW in the fornA 

(e£ + + e£(z) + W"(zty tf"(p,*) = E^(p,z), (2) 



where the isotropic kinetic energy is given by 

! _|_ „2 

(3) 



Px+Py 



P 2m, ' 

which includes the effective mass m,. The £%(z) energy is £™(z) = U(z) in the wells, and 
£ c (z) = AE C + U(z) in the barriers, where AE C is the band offset for c conduction band. 
Whenever energy values are less than AE C underbarrier penetration (and tunneling) will 
be permitted and described by the boundary conditions. Now, parabolic approximation is 
justified because energy values under consideration are smaller than the gap energy e g in 
the narrow region. Band diagram for ACQW is shown in Fig. 1. The potential U(z), for an 
uniform transverse electric field, is U (z) ~ eF±z. 

The magnetic energy W ll (z) 1 for not very strong magnetic fields, is described by 

W"{z) = -V>{z) [a x P ] 2 + ^-fiBH&y, (4) 

where \xb = \e\h/(m e c) is the Bohr magneton and a is the Pauli matrix. Finally, the 
characteristic spin velocity V^(z), and the effective (7-factor are 



V»{z) = J_^(z)/dz t _ m 



4m M e g 2m, 



1 + z cV " — 



£ 9 



(5) 
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FIG. 1: (Online color) Conduction band diagram for undoped asymmetric double quantum well. 
Horizontal thin lines show electron energy levels and thin curves correspond to squared wave 
functions close to the resonance. 

The potential of interfaces determines a part of the spin dependent contribution through 
the parameter \. Actually, % takes into account the spin-orbit coupling due to the abrupt 
potential of the heteroj unction at each interface, as2*Ai 

2 ^ 



X = - J m^{z) dz, (6) 

where the integral is taken over the width of the abrupt interface 25. 

In the present case, where U(z) is linear with z, we get de£(z)/dz = dU(z)/dz = eF±, 
which does not depend on z. Thus, the external applied electric field F± shapes the spin 
velocity V^{z) = and the (/-factor, g^{z) = We can take the characteristic spin 

velocity for each layer as 

4m (1 e 3 ' ^ 
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and the abrupt interface parameter as 

2eF ± 5 + AE C 



X = 



2e r , 



AE C 

2e r . 



(8) 



Lastly, we introduce the Pauli splitting energy wh, caused by the magnetic field, as = 
(g»/2)pi B H. Thus, Eq. (4) becomes 

W» = v»laxp] z + w» H a y , (9) 

which has lost the z dependence. 

Because and e£ do not depend on z, we can factorize fundamental solutions of Eq. 
(2), \l> M (p,z), as products of ^^(p) an d (p^ a (z) functions, with a = ±1. The a value refers 
to the two possible spin orientations. For an ACQW under a transverse electric field F±, 
the p-dependent spinors tp^ip) can be obtained from 



(10) 



in the form 



^(p) - T2 



(v fl p + + w^)/iw£ 
(v>_ +w II H )/iw£ 



where 



(V"p x + O 2 + (V%) 2 , 



(11) 



(12) 



and the energy quasi-paraboloids are 

e^ = e^ + (7w^, (13) 

Now, z-dependent functions <f^ a (z) are obtained from the second order differential equa- 
tion 



2m„ 



(14) 



For an ACQW under a transversal electric field F±, eigenstate functions of Eq. (14) are the 
known linear combination of the Airy Ai- and .Bi-functions. Thus, the general solution of 
Eq. (2) can be written as 



^ + (p,z) = ^ + (p) [a, + A^{z) + , 

^-(p, Z ) = r~(p) [a,-Ar(z)+b^Br(z)} , 



(15) 
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where a^, b, a are four constants by region to be solved through the interface conditions of 
continuity of the wave functions and current. Because there are two wells, and three barriers 
(Fig. 1), there are four interfaces and four boundary conditions by interface. 

Next, to simplify the calculation of the wave functions and the dispersion relations of the 
electronic levels through the boundary conditions*^, we create two auxiliary parameters: a 
length and an energy 

% 2 \ 1/3 „ ^ 



and a set of momentum dependent functions 

Pp— ■ m 

IWp 
P+V M + Wjj 



pp+ 



■ A 1 
ZWp 



, M+ V" (Pi + Pi) + WhP* 

Zip X" 



where 



rn+ _ w hPv 

/ 2 M P " = -/r P + (17) 



P+ = Px + ipy, and p- = p x - ip y . (18) 



Next, using preliminary quasi parabolic dispersion relations we construct the Airy 
function arguments 

where 5^ acts as a Kronecker function: = 1 when fi = b, and 5£ = when fi = w. 

We use the transfer matrix method to obtain wave functions and energy dispersion rela- 
tions of the electronic levels. In the present case, we have improved the standard method 
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by using 4x4 matrices at each interface, whose elements are 



MZ 2 
M1 3 
Ml, 
M» 21 
M» 2 

M^ 

Ml 
Ml 
Ml 
M^ 



E,p 
E,p 
E,p 
E,p 
E,p 
E,p 
E,p 
E,p 
E,p 

E,p 

E,p 

E,p 

E,p 

E,p 

E,p 

E,p 



Ai (e P+ ) , 
Me P -), 

^m&)+fgf&.Bi (e P+ )), 

[p" P -Bi'(^-) + f£Bi{e P -)] , 
[p" P+ Aif(e P+ ) + ft* (e P+ )} , 
[p" P+ Bi'(e P+ )+f^Bi (&.)], 
[Ai'i^+f^P^-Ai (e P -)}, 
[Bi'(e P -) + tt pP » P -Bi\ (e P -) , 



rn 

m 
rn 

m 
m 

rn 
rn 

rn 
m 

m i 
rn 

m i 
m 

rn 
m 



m 



(20) 



where Ai'(^ p ) means dAi(££ p ) / dz and Bi'(^ p ) means dBi(^ p )/dz. Now we are ready to 
generate transfer matrices, S^(z : E, p) = M^(z, E, p). Finally, electronic levels for each 2D 
momentum p = (p x ,p y ) are obtained from S44 (E, p) = 0, where 



n -1 



S (E, p) = ]S b (L , E, p)j • S W (L , E, p) • [^(Lx, E, p) 
S\L U E, p) • \S\L 2 , E, p)l ~* • S W (L 2 , E, p) • 



S w (L 3 ,E,p) ■ S\L 3 ,E,p 



(21) 



In the above matrix product, z = Li denotes interfaces position in the growth direction, 
starting from the left side. 
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Calculations give us the two coupled spin up paraboloids E k+ (p) and spin down 
paraboloids E k _(p), where k = 1,2 corresponds to the deepest levels of the ACQW. Once 
obtained the involved coefficients a^, we normalize wave functions for each momentum 
P 

The scheme of Fig. 1 includes the two resonant energy levels and the respective wave 
functions for p = 0. Although there are four levels only two are observable in this figure. 
This is because spin sublevel splitting is much smaller than electronic level energy distance 
and differences between spin down and spin up wave functions are not visible at p = 0. 

Finally, the density of states can be obtained by using the well-known expression 

^) = E/(2^p 5 ( £ -^(P))- ( 22 ) 

Peculiarities of p(e) can be analyzed through the photoluminescence excitation (PLE) in- 
tensity for the case of near-edge absorption, I ple, because both quantities are related hy-<^ 

I ple ~ ^2 |e ■ v cv \ 2 5(e Xc - £\ v - hu) ~ p(hAuj) (23) 
for very low temperature or 

I ple ~ ' e ' v ™\ 2 G(e Xc - e Xv - hw), (24) 



when one includes electron-phonon scattering. In the above expression, the Gaussian func- 
tion is G(x,xq) = ^^exp — (^^j ■ The Gaussian halfwidth 7 is related to the scat- 
tering and relaxation processes and, thus, to the temperature- 20 . Expressions (23, 24) are 
valid provided the interband velocity v c „ does not depend on in-plane momentum. Here e 
is the light polarization vector, Aw = uj — e g /h, and e Xc , e Xv are the conduction and valence 
band levels, respectively. 



III. RESULTS 



Let's start this section with numerical results for In x Ga\- x As — In y Ali_ y As-ba,sed AC- 
QWs, with x = 0.53 and y = 0.52. We have chosen this particular structure because we 
have reliable data for basic parameters^. We have considered two InGaAs wells of 70 and 
100 A wide separated by a 20 A InAlAs barrier. We have also applied an electric field of 30 
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kV/cm, which corresponds to a spin velocity v™ = 2.6 x 10 5 cm/s for the InGaAs QWs, and 
v 6 = 1.4 x 10 5 cm/s for the InALAs barriers, with a transition spin velocity region across 
de abrupt interface. This electric field is slightly higher than needed to achieve resonance 
between the deepest levels of both wells (28 kV/cm). To calculate interface contributions 
we have used a typical abrupt interface size of 5 ~ 3A for InAlAs — InGaAa^-. We have 
also applied in-plane magnetic field of 0.01 T. 

Figures 2 (a — d) show normalized squared wave function for p y = 0, versus z and the 
dimensionless momentum p x /po, where po = m w \i w . Upper panels (a, h) correspond to the 
first deepest level for the two different spin orientations. As expected for an electric field 
beyond the resonance, charge density is mainly located in the left narrow QW. Consequently, 
the opposed happens for the second resonant level as can be seen in the lower panels (c, d) 
Analyzing the behavior of wave functions versus momentum p x and spin orientations by 
comparing panels (a) and (b), a new peculiarity appears. While the charge distribution 
coincide for both down and up spins at the zone center (p x = 0), there is a charge transfer 
between wells for increasing \p x \. For spin down case [panel (a)] charge experiences a little 
shift from left narrow well to the right wide one. As a result, to conserve the charge, the 
opposite occurs for spin up electrons [panel (b)} where charge goes from the right to the 
left well. For the higher resonant level we find a misleading ambiguity. It might seem the 
behavior is the opposite to the previous one because now, the spin down case [panel (c)] 
shows a charge transfer from right to left QW as \p x \ increases and, conversely, for spin up 
electrons [panel (d)}. However, considering the relative charge concentration between wells, 
we can realize there is a similar behavior for both resonant levels. For spin down electrons 
there is a charge shift from well with higher concentration to the other well [panels (a), (c)] 
and conversely for the spin up electrons [panels (6), (d)]. The reason for this charge shift 
lies in the energy term w£, which induces a breaking of the momentum symmetry. Because 
of w£ is an essential part of the argument of the Airy functions, the behavior of the wave 
functions is significantly affected. 

Fig. 3 shows the near parabolic dispersion relations of the two coupled levels and their 
corresponding spin up and down sublevels, for the electric and magnetic fields under con- 
sideration. It can be seen the p y paraboloids symmetry agrees with Eqs. (11) and (18). 
Although there is a little difference between spin paraboloids, due to the large energy dif- 
ference between the resonant levels of both wells (~12 meV) and the splitting of the spin 
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FIG. 2: (Online color) Squared wave function vs normalized p x /po for the first spin resonant levels. 

sublevels (~0.01 meV at p y = 0) it is not possible to distinguish minima behavior. Thus, we 
have enlarged in Fig. 4 the bottom of the pair of paraboloids (spin up and spin down) for 
the ground level. After resonance this level would be associated to the deepest level of the 
left narrow QW, if wells were decoupled. As expected, both paraboloids shift in opposite 
p x directions resulting in a sublevel anticrossing. 

Considering a p y = section of the former figure we can get a more accurate 2D repre- 
sentation (Fig. 5) of the anticrossing, minima p x position, and energy splitting. The inset 
displays anticrossing area enlargement. In order to have a more detailed overview of the 
anticrossing region we have also included in Fig. 6 the contour plot around anticrossing for 
different constant energy values. 

Next, we analyze the density of states p (e) because this function is proportional to 
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FIG. 3: (Online color). Near parabolic dispersion relations of the two coupled levels close to the 
resonance, and their corresponding spin up and spin down sublevels. 

the photoluminescence excitation (PLE) intensity—. And PLE is one of the most used 
techniques to get information of quantum structures. In this case we have used magnetic 
field intensities varying from to 0.1 T. The shape of p (e), shown in Fig. 7, is characterized 
for the peaks coming from he e~ 1//2 -divergence corresponding to the zero-field case. As 
expected, these peaks disappear gradually by growing magnetic field. In turn, interfaces 
contribute with a manifested delay in the quenching of the p (e) singularities, as well as an 
additional broadening of these peaks. That is, although the singularities should only appear 
at zero magnetic field, they still remain at the band anticrossing position for low magnetic 
fields, as shown in Fig. 7. Another significant feature is that the p (e) anticrossing peak is 
softened when increasing barrier height, disappearing for lower magnetic fields than used in 
this work 3 

Finally, Fig. 8 shows PLE spectra for different Gaussian halfwidth 7 at a fixed H = 0.1 
T. For the magnetic field used before (H = 0.01 T) the two adjacent peaks, corresponding to 
each resonant pair of states, overlap. Thus, we have used a magnetic field ten times higher 
because this field allows us to tell the peaks apart for small 7 values. Evolution of the first 
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FIG. 4: (Online color) Bottom of the two deepest paraboloids just after resonance and detail of 
spin up and spin down anticrossing. Levels are mainly located in the narrow (left) QW. 

two peaks with 7 is depicted in Fig. 9. As can be seen, PLE peaks corresponding to the 
two different spin transitions are still distinguishable for 7 values beyond 1 meV. 

IV. CONCLUSIONS 

In this paper we have analyzed the electron spin behavior in narrow-gap ACQWs under 
transverse electric and in-plane magnetic fields. We used the Kane model with nonsymmet- 
ric boundary conditions, caused by the two different wells width, to solve the eigenvalue 
problem. Based in this model and the transfer matrix approach we have performed a use- 
ful tool to tackle any layered structure with abrupt interfaces and subjected to different 
perturbations. The model allows us the study of the spin peculiarities of level anticrossings. 

Because interface contributions oppose intrinsic spin-orbit effect, mechanisms that mix 
the Pauli contribution with the two kinds of spin-orbit contributions (from a low magnetic 
field and from heterojunctions) are different. As a result numerical calculations lead to 
magnetoinduced variations of the energy spectra under in-plane magnetic field strength 
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FIG. 5: (Online color) Two dimensional dispersion relations for p y /po = 0. Inset shows an magnified 
image of the spin up and spin down anticrossing. 

even lower than Tesla for typical parameters of InGaAs/InAlAs structures. 

There are essentially two techniques to check electron energy spectrum changes obtained 
in the present work. One of them, the magnetotransport measurements of Shubnikov- 
de Haas oscillations under nearly in-plane magnetic fields- 13 -'— provide us information of 
in-plane magnetoinduced modifications of electron energy. The other technique to find 
characteristics of energy spectrum is the mid-infrared PLE spectroscopy—. This technique 
is suitable because it provides direct information of the energy spectrum when interband 
transitions are modified by in-plane magnetic fields. Moreover, PLE intensity is directly 
related to the density of states, which we can evaluate with the present model from the 
energy dispersion relations. 

A detailed analysis beyond the Kane model does not modify essentially our formulation 
of the eigenstate problem. Thus, the main conclusions of the paper on the changes of slow- 
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FIG. 6: (Online color). Contour plot of the deepest energy levels around anticrossing. (a) E\„ = 
83.6100 meV; (b) E la = 83.611 meV; (c) E Xa = 83.612 meV; and (d) E l(J = 83.614 meV. 

electric field and hetero junction-induced spin-splitting mechanisms under in-plane magnetic 
field for electron spin tunneling between wells in ACQWs, remain still valid. We have also 
used along the work the assumption that the confined potential is well described. Thus, the 
effective transverse field F± provides correct estimations both for magnetoinduced changes 
and for the character of the dispersion laws. However, more detailed numerical calculations 
are needed to describe the kinetic behavior mentioned above. We will return to this point 
in a forthcoming work where we will analyze the spin dynamics in ACQWs. In summary, 
we have shown that essential variations of the dispersion relations in ACQWs occurs when 
we include contribution of the abrupt barriers. Electron density of states and PLE reflect 
these modifications in narrow-gap structures under in-plane magnetic fields when tunneling 
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FIG. 7: (Online color) Density of states for the same electric field of previous figures and different 
magnetic fields. 

between wells occurs. 
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